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1. Introduction 

The thermodynamics of strongly interacting matter at finite temperature has been studied ex- 
tensively. It has been shown that the QCD transition is a cross-over [jlj]. This result is of importance 
for the physics of the early universe. It is also interesting to examine the QCD phase diagram at 
finite chemical potential \i% ^ which would be relevant for the description of heavy ion colli- 
sions. Lattice simulations in this regime have been hindered by the sign problem ever since the 
first attempts [gj. The sign problem originates from the fact that the fermion matrix M loses its 
75 hermiticity at finite chemical potential M ^ 75^75. Thus the fermion determinant det(M) be- 
comes complex and spoils any method based on importance sampling. There are several different 
methods to circumvent this problem for small }Xb based on either performing simulations at imag- 
inary jig and then analytically continuing the results to the physically interesting case of real }i. B 
[^] or reweighting samples produced with a real and positive det(M) [Q, ||]. For our calculations 
we used a reweighting technique, where all observables measured on the = configurations are 
expanded in }Xb- The analysis was performed using the same configurations as in [0]. 



2. Reweighting and Taylor expansion 

The reweighting technique uses configurations generated at /ig = (where the fermion deter- 
minant is real and positive), to access observables G at finite chemical potential. 
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The expectation value (G) at /ig and j8 is rewritten as the ratio of two expectation values at /ig = 
and j8o- The chemical potential range where the reweighting is effective is limited by the decreasing 
overlap between the sample generated at /ig = and the target ensemble at )J.b 7^ 0. As the distri- 
bution of the sample becomes narrower with increasing volume the overlap problem gets also more 
severe, thus requiring higher statistics when approaching the continuum limit. It is quite expensive 
to calculate the determinant det(M(/ig)) exactly, especially for large lattices. For that reason we 
decided to use a Taylor expansion in ju. This technique has been used several times, especially by 
the Bielefeld-Swansea Collaboration (for a recent paper with sixth order expansion see [f7|]). For 
gluonic observables which do not depend on the chemical potential explicitly this expansion is very 
easy to perform. 
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One can also note, that at ]X q = the odd order derivatives of lndet {M{}X q )) are purely imaginary 
and the even order derivatives are real. This follows from the fermion matrix identity: 

M^{pL q ) = y 5 M{-pL q )y 5 (2.5) 

Therefore it is clear that for a real observable the linear term in the expansion is zero. 

3. Observables 

All observables that we compute in this work namely the Polyakov loop and the strange quark 
number susceptibility do not depend explicitly on }Xb and can thus be expanded like gluonic ob- 
servables. 

3.1 Polyakov loop 

The Polyakov loop is given by the trace over a product of gauge links along a line in time 
direction: 

P = 7^3 L tr [^4 (x, 0) 04 (x, 1) . . . U 4 (x,N t - 1)] (3. 1) 

In pure gauge theory a change in the expectation value of the Polyakov loop to a non-vanishing 
value signals the spontaneous breakdown of the Z(3) symmetry and thus the onset of deconfine- 
ment. The Polyakov loop is also related to the quark-antiquark free energy at infinite separation: 

|(P)| 2 = exp(-A^(r^oo)) (3.2) 

Here AF qq denotes the difference of the free energies of the quark-gluon plasma with and without 
the quark-antiquark pair. It is then possible to renormalize the Polyakov loop by renormalizing the 
free energy of the quark-antiquark pair. For that purpose one uses the renormalization condition 
Vii( r o) = to renormalize the potential at T = for each lattice spacing As the ultraviolet 
divergencies are the same at T = as at finite T the same renormalization shift can then be used to 
renormalize the free energy. The renormalized Polyakov loop is given by 

\(P R )\ = KF)|exp(y(r )/(2r)), (3.3) 

where V(ro) is the unrenormalized T = potential obtained from Wilson-loops. The transition 
temperature can be defined as the peak in the temperature derivative of the Polyakov loop, that is 
the inflection point of the Polyakov loop curve. 

We also note that the linear term in the Taylor expansion of the Polyakov loop is zero. This 
can be shown by using the CPT symmetry It is enough to show that (P)^ = (P)!. This follows 
from the CPT transformation of the expectation value of the Polyakov loop: 

(P)^ = (p)™ = J W CPT e- s ^ f/CPT )det(M M ([/ CPT ))P CPT 

= J &Ue~ s ^ [det(M M (£/))] * P* = ((*%)* (3.4) 
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3.2 Strange quark number susceptibility 

The strange quark number susceptibility % s is defined by 

% s 1 d 2 \ogZ 



T 2 TV dpi} 



(3.5) 

n s =o 



It can be related to the event-by-event fluctuations in heavy-ion experiments. As % s has a well 
defined continuum limit, renormalization is not necessary. The transition temperature can be de- 
fined as the peak in the temperature derivative of the strange quark number susceptibility, that is 
the inflection point of the susceptibility curve. 

As the strange quark number susceptibility is a real observable the linear term in the \Ir ex- 
pansion vanishes as was shown above. 



4. Simulation Setup 

We used a Symanzik improved gauge and stout-link improved staggered fermionic lattice ac- 
tion in order to reduce taste violation The configurations were generated with an exact RHMC 
algorithm. We determined a line of constant physics (LCP) using physical masses for the light 
quarks m u ^ as well as for the strange quark m s . The LCP was fixed by setting m^/ '/k and mK/m K 
to their physical values We used four different lattice spacings N t = 4,6,8, 10 with an aspect 
ratio of N s /N t = 4. The scale was fixed by fx and its unambiguity checked by calculating rriK*, 
f n and ro. The operators 0\ and were computed using a random noise estimator for the traces. 
The number of random vectors was chosen to give an error of the same size as the statistical error. 
It was between 50 and 160 in all cases. 



5. Results 
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Figure 1: On the left hand side the temperature dependence of the expectation value of the Polyakov loop 
is shown for fia = Q and ju a = 0. 1 . The graph on the right hand side depicts the temperature dependence of 
the strange quark number susceptibility for jj. a = and jia = 0. 1 . The results were obtained on the N t = 6 
lattices. 
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Figure 1 shows the Polyakov loop and the strange quark number susceptibility for \i a = and 
/I a = 0.1 at N t = 6. We want to compute the curvature K of the phase transition line at jj. B = 0, 
which we define by 



T c (Hb 



that is 
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(5.2) 



H B =0 



Here it is assumed that the phase transition line can be well described by a purely quadratic term in 
{JL B , whereas the fourth order term is negligible. Now we use 
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to directly extract the curvature from d 0/dfig and d£?/dT. Here G stands for the Polyakov loop or 
the strange quark number susceptibility respectively. Figure 2 shows —T c dT c {\lB) I d\l B computed 
by dT c (pL B ) i 'dflg = d% s /dil B l (d% s /dT)~ for the N t =41attices. We use an average over the points 
near T c {}JLb = 0) to extract the curvature at T c and can thus increase the statistics. 
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Figure 2: This plot shows the -T c dT c {jX B )/djxl computed by dT c (/j. B )/dfi]j = dx./d^idXs/dTy 1 for 
the N t =4 lattices. 

In our case d@ jd\l\ is directly given by the quadratic term in the \i expansion. Whereas 
d<£> '/dT can be extracted from a fit to the Polyakov loop and the quark number susceptibility near 
T c . For one point we also calculated d6 /dpL B by exact determinant evaluation JTo| ] and found a 
perfect agreement. 

We determine the error of the curvature by carrying out the whole analysis for different 
jackknife-samples and by varying the fit range. 

To extract the continuum limit of K one needs at least three points in the scaling region. It was 
established in ^ that our N t =4 lattices aren't in this region. Thus we use the N t = 6,8,10 results 
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Figure 3: The left plot shows the phasediagram in the pig-T plane for the lattice spacing N t = 4. The red 
curves give the crossover line computed by finding the inflection point of the Polyakov loop, the black one 
corresponds to the computations using the strange quark number susceptibility. The big errorbars at the start 
of the curve denote the total error of T c (ji = 0), the smaller errorbars include just the error of the shift AT. 
As the phasediagram was extracted from a Taylor expansion it is only correct to leading order. On the right 
hand side the lattice spacing dependence of the curvature K is displayed. No safe continuum extrapolation is 
possible with our present statistics. Here the red and black lines correspond to the results obtained from the 
Polyakov loop and the strange quark number susceptibility respectively. 



for the extrapolation, see Figure 3. Unfortunately the present statistics for the N t = 10 lattices is 
very limited and thus the errors are large. It is an ongoing project to increase the statistics to give a 
reliable continuum value for the curvature. 



6. Summary 

We have computed the curvature K of the pL-T pseudocritical line at /j.b = 0. The simulations 
were performed at physical quark masses with nf = 2+1 flavours for 4 different lattice spacings. 
We used a Taylor expansion in /ig of the Polyakov loop and the strange quark mass susceptibility 
and determined their inflection points in order to locate the crossover line. 



Acknowledgments 

The computations were carried out on the 370 processor PC cluster of Eotvos University, on 
the 1024 processor PC cluster of Wuppertal University, on the 107 node PC cluster equipped with 
Graphical Processing Units at Wuppertal University JTT| ] and on the Blue-Gene/L at FZ Julich. We 
used a modified version of the publicly available MILC code [|12|] with next-neighbour communi- 



cation architecture for PC-clusters [13]. Partial support of grants of DFG F0 502/1, EU I3HP and 
OTKA AT049652 is acknowledged. 



6 



Curvature of the phase transition line in the pt-T plane 



Christa Guse 



References 

[1] Y. Aoki, G. Endrodi, Z. Fodor, S. D. Katz, K. K. Szabo, The order of the quantum chromodynamics 
transition predicted by the standard model of particle physics, Nature 443 (2006) 675-678 

[hep-lat/0611014] 

[2] I. Barbour et. al., Problems with Finite Density Simulations of Lattice QCD, Nucl. Phys. B275 (1986) 
296 

[3] Ph. de Forcrand, O. Philipsen, The QCD phase diagram for three degenerate flavors and small baryon 
density Nucl. Phys. B673 (2003) 170-186 [hep-lat/0307020] 

[4] I. M. Barbour, S. E. Morrison, E. G. Klepfish, J. B. Kogut, M. P. Lombardo, Results on finite density 
QCD Nucl. Phys. Proc. Suppl. 60A (1998) 220-234 [hep-lat/ 97 05042] 

[5] Z. Fodor, S. D. Katz, A new method to study lattice QCD at finite temperature and chemical potential, 
Phys. Lett. B534 (2002) 87-92 [hep-lat/0104 001] 

[6] Y. Aoki, Z. Fodor, S. D. Katz, K. K. Szabo, The QCD transition temperature: Results with physical 
masses in the continuum limit, Phys. Lett. B643 (2006) 46-54 [hep- 1 at/0609068] 

[7] C. R. Allton et. al., Thermodynamics of two flavor QCD to sixth order in quark chemical potential, 
Phys. Rev. D71 (2005) 054508 [hep-lat/0501030] 

[8] B. Alles, E. M. Moroni, A proposal for simulating QCD at finite chemical potential on the lattice, 
(2002) [hep-lat/0206028] 

[9] Y. Aoki, Z. Fodor, S. D. Katz, K. K. Szabo, The equation of state in lattice QCD: With physical quark 
masses towards the continuum limit, JHEP 01 (2006) 089 [hep-lat/0510084] 

[10] Z. Fodor, S. D. Katz, Lattice determination of the critical point of QCD at finite T and mu, JHEP 03 
(2002)014 [hep-lat/0106002] 

[11] G. I. Egri et. al., Lattice QCD as a video game, Comput. Phys. Commun. 177 (2007) 631-639 

[hep-lat/0611022] 

[12] MILC Collaboration, public lattice gauge theory code, see http://physics.indiana.edu/sg/milc.html 

[13] Z. Fodor, S. D. Katz, G. Papp, Better than 1 jMflops sustained: A scalable PC-based parallel 

computer for lattice QCD, Comput. Phys. Commun. 152(2003) 121-134 [hep-lat/0202030] 



7 



